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ABSTRACT 

A new, three-dimensional, shock capturing, hydrodynamic code is utihzed 
to determine the distribution of hot gas in a CDM+A model universe. Peri- 
odic boundary conditions are assumed: a box with size 85/i~^Mpc, having cell 
size 0.31/i~-^Mpc, is followed in a simulation with 270^ = 10^-^ cells. We adopt 
n = 0.45, A = 0.55, h = i^/lOOkm/s/Mpc- 0.6, and then, from COBE and light 
element nucleosynthesis, erg = 0.77, Clb = 0.043. We identify the X-ray emit- 
ting clusters in the simulation box, compute the luminosity function at several 
wavelength bands, the temperature function and estimated sizes, as well as the 
evolution of these quantities with redshift. This open model succeeds in match- 
ing local observations of clusters in contrast to the standard Q — 1, CDM model, 
which fails. It predicts an order of magnitude decline in the number density of 
bright {hu = 2 — lOkeV) clusters from ^ = to ^ = 2 in contrast to a slight 
increase in the number density for standard = 1, CDM model. 

This COBE-normalized CDM+A model produces approximately the same 
number of X-ray clusters having > lO^^erg/s as observed. The background 
radiation field at IkeV due to clusters is approximately 10% of the observed 
background which, after correction for numerical effects, again indicates that the 
model is consistent with observations. 

The number density of bright clusters increases to 2; ~ 0.2 — 0.5 and declines, 
but the luminosity per typical cluster decreases monotonically with redshift, with 
the result that the number density of bright clusters shows a broad peak near 
z = 0.5, and then a rapid decline as 2; — > 3. The most interesting point which 
we find is that the temperatures of clusters in this model freeze out at later 
times {z < 0.3), while previously we found in the CDM model that there was 
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a steep increase during the same interval of redshift. Equivalently, we find that 
L* of the Schechter fits of cluster luminosity functions peaks near z — 0.3 in 
this model, while in the CDM model it is a monotonically decreasing function 
of redshift. Both trends should be detectable even with a relatively "soft" X-ray 
instrument such as ROSAT, providing a powerful discriminant between Q = 1 
and n < 1 models. Detailed computations of the luminosity functions in the 
range Lx = 10^*^ — lO^^erg/s in various energy bands are presented for both 
cluster cores (r < 0.5/i~-'^Mpc) and total luminosities (r < lh~^Mpc). These are 
to be used for comparison with ROSAT and other observational data sets. They 
show the above noted negative evolution for the open model. 

We find little dependence of core radius on cluster luminosity and the de- 
pendence of temperature on luminosity log kT^ = A + B log Lx, which is slightly 
steeper {B — 0.32 ± 0.01) than indicated by observations {B — 0.265 ± 0.035), 
but within observational errors. In contrast, the standard = 1 model predicted 
temperatures which were significantly too high. The mean luminosity- weighted 
temperature is l.SkeV, dramatically lower (by a factor of 3.5) than that found 
in the Q = 1 model, and the evolution far slower (—30% vs —50%) than in the 
Q = 1 model to redshift z — 0.5. A modest average temperature gradient is 
found with temperatures dropping to 90% of central values at OAh~^Mpc and 
to 60% of central values at 0.9/i~^Mpc. 

Examining the ratio of gas-to-total mass in the clusters, we find a slight an- 
tibias (b — 0.9 or {jiff)ci — 0.083 ± 0.007), which is consistent with observations 
[(^^)o6s = 0.097 ± 0.019 for the Coma cluster for the given value of h, cf. , 
White 1991]. 

Cosmology: large-scale structure of Universe - hydrodynamics - Radiation Mech- 
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anisms: Bremsstrahlung - X-ray: general 
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1. INTRODUCTION 

In the preceding two papers of this series (Kang et al. 1994 = KG OR 
hereafter; Bryan et al. 1994), we have examined the properties of X-ray clus- 
ters in the standard COBE- normalized CDM universe and reached two primary 
conclusions: 1) the standard CDM model overproduces bright X-ray clusters 
{Lx > lO'^'^erg/s) by a factor in excess of five; 2) observations of the ratio of 
gas-to-total mass in great clusters of galaxies, combined with our simulations, 
impliy that we live in an open {fl = 0.2 to 0.3) universe (with or without a 
cosmological constant). There are other well-known difficulties with the stan- 
dard CDM scenario (c/. Ostriker 1993 for a review), and there are advantages to 
open CDM models (c/. Efstathiou 1992) which arise from other considerations, 
independent of cluster X-ray properties. Motivated by this knowledge, we turn 
in this paper to examine an open CDM model with a cosmological constant. 
We wish to investigate the difference in the evolutionary behavior as well as the 
z = properties of X-ray clusters in these two different models. There is every 
expectation that cluster X-ray properties should provide a strong discriminant 
among cosmic theories. We will, throughout this paper, make comparisons with 
KCOR wherever it is possible, since that calculation was made with identical 
physical assumptions and with identical numerical modeling techniques adopted. 

In §2 we outline the method and initial conditions, in §3 the results, and in 
§4 we assemble our conclusions. 

2. METHOD AND INITIAL CONDITIONS 

2.1 Method 

The superiority of the new TVD code over conventional hydrodynamic codes 
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(e.g., Cen 1992) and various tests of the code have been presented in Ryu et 
al. (1993); essentially the shock capturing technique improves resolution by 
approximately a factor of 2 — 3 for a given grid (i.e., nominal resolution). 

The simulation reported on in this paper did not include any atomic pro- 
cesses, i.e., no cooling or heating was added, except for the adiabatic cooling 
due to the general expansion of the universe, and "heating" occurs only due to 
adiabatic compression or to entropy generation at shock fronts. For the hot gas, 
which we will discuss in this paper, this approximation is valid, since the cooling 
time exceeds the Hubble time by a fair margin. One can, after the fact, compute 
the emissivity with allowance for line emission. Doing this would increase the 
computed luminosity significantly at and below IkeV but would have little effect 
in the 2 — lOkeV band. We will return to this important matter in a subsequent 
paper devoted to a comparison with observations. Here we wish, primarily, to 
contrast open, f2 < 1, and standard = 1, variants of the CDM scenario. 

2.2 Initial Conditions 

We adopt a CDM+A model with the following parameters: n = 1, /i = 0.6, 
Q, = 0.45, Qb = 0.043, A = 0.55, and ag = 0.77. Note that the amplitude normal- 
ization of the power spectrum is determined by COBE observations (Efstathiou, 
Bond, & White 1993; Kofman, Gnedin, & Bahcall 1993) parameterized by ag 
to translate into conventional notation. This chosen value of as corresponds to 
a "bias" of 6 = l/ug = 1.30, close to what is physically expected on this scale 
(c/. Cen & Ostriker 1992,1993) lending further credibility to the adopted model. 
Our box size is 85/i~^Mpc with N — 270^ cells and 135"^ dark matter particles, 
so our nominal resolution is 0.31/i~^Mpc with our real spatial resolution slightly 
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worse than this. The box size is determined by a compromise between two con- 
siderations. A larger box would allow longer waves and give us more of the rare, 
high temperature, high luminosity clusters. A smaller box (with a fixed N) would 
allow us to resolve the cluster structure better. The choice of Clb is consistent 
with light element nucleosynthesis (Walker et al. 1991). The power spectrum 
transfer function is computed using the method described in Cen, Gnedin, & 
Ostriker (1993). Gaussian initial conditions are used and the same set of random 
numbers adopted as in KCOR. 

3. RESULTS 

The X-ray clusters in the simulation are identified as follows. We first 
calculate the total X-ray luminosity due to thermal Bremsstrahlung (assuming 
primeval composition and neglecting lines) for each cell, given the cell density 
and temperature, assuming that hydrogen and helium are fully ionized (which is 
always true in the regions like great clusters of galaxies). The detailed formulae 
were presented in KCOR. Note that in the following discussion all units of length 
are given in co-moving, not metric, coordinates. 

First we tag all cells having total X-ray luminosity higher than lO^^erg s~^. 
This emissivity corresponds to 3.2 x 10^^erg/s//i~'^Mpc'^, which is 3.2 times the 
mean box emissivity at z=0, and is 9.1 times the mean at z=3. The number 
of X-ray bright cells defined in this way is 47,826, which comprises a fraction 
2.4 X 10~^ of the box volume. These are selected as X-ray bright cells. Then we 
find the local maxima (by comparing Lff of each X-ray bright cell with that of 26 
neighboring cells) and identify them as the centers of the X-ray clusters. Having 
defined the centers of the X-ray clusters, we go back to the whole simulation box 
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to define our X-ray clusters. We analyze the simulation in the following two dif- 
ferent ways (which correspond to spheres of radius 0.5h~^Mpc and 1.0/i~^Mpc). 
First, each cluster core consists of 27 cells (26 cells surrounding the central cell 
plus the central cell). These 27 cells are weighted so that the total volume of the 
cluster equals the volume of a sphere of radius 0.5h~^Mpc as appropriate for ob- 
servationally defined X-ray clusters (see KCOR for details). A similar algorithm 
is used for the 1.0/i~^Mpc volumes. 

Our (85/i~^Mpc)^ box a,t z — contained (0,0) clusters with (total, core) 
luminosity brighter than lO^^erg/s, (1,0) brighter than lO^^erg/s, (10,5) brighter 
than lO^^erg/s, (66,40) brighter than lO^^erg/s and (257,174) brighter than 
lO'^^erg/s. We explain the absence of clusters brighter than lO'^^erg/s as sim- 
ply due to the size of our box. An additional factor of at least two in scale (and 8 
in computer resources) would be required to significantly improve on the quoted 
results. An immediate comparison to the KCOR results is possible if one notes 
that for Lx > lO'^^erg/s, the KCOR paper found (8,1) clusters with (total, core) 
luminosity, higher than the specified limit compared to (1,0) in this work. The 
difference is significant at this level of luminosity and we will see below (Figure 
lb and 4b) that, while the standard CDM model overproduces bright X-ray clus- 
ters, this open model provides an adequate fit to observations. The fraction of 
the box in brightest cells which provide (50%, 90%) of the total box emissivity is 
(3.8 X 10-^ 8.8 X 10-^) 

It is convenient to fit the luminosity function to the three parameter Schechter| 
function 



n{L)dL = no{L/L*ye-^/^'d{L/L*) 
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(1) 



Luminosity functions for cluster cores were computed in four frequency bands: 
total (bolometric) luminosity, 0.3 — 3.5keV, 0.5 — 4.5keV, and 2 — lOkeV, and are 
displayed in Figures la, 2a, 3a and 4a, respectively. The results for entire clusters 
(emission from a lh~^Mpc sphere) are presented for the same frequency bands 
in Figures lb, 2b, 3b and 4b. The figures show the domain of cluster properties 
in which the observations and our computations overlap most: 10^°erg/s< < 
lO'^^erg/s and < z < 1. Also shown (cross-shaded areas) in Figures (lb) and 
(4b) are the observations from Henry & Arnaud (1991) and from Henry (1992), 
respectively. The comparison shown in (4b) is the more reliable one, since line 
processes (omitted in our computation) are unimportant in the 2-lOkeV band. 

We see that the computed number densities of bright clusters {h'^L > 10^'^erg/s 
are in accord with the observed ones while in KCOR we found that those of the 
CDM model are above the observed mean by a factor of about 5. We have 
computed approximate Schechter function fits to the results, with the numerical 
parameters {no,L*,a) as a function of redshift collected in Table (1), and the 
simulated data extended to z = 5. Also in Table (1) we integrate over the cluster 
luminosity function, using the Schechter fit, jd = noI'*r(2 — a), showing the 
results in the second-to-last column. We give also, in the last column, the total 
emissivity from the box as jgas, which includes the emission from lower density 
regions further from cluster cores than 1.0/i~^Mpc and also from clusters whose 
central emissivity is less than our cutoff value. Note that L* and L are in units 
of lO^^erg/s; no, n(L > lO^^) and n{L > 10^^) are in units of lO'^/i^Mpc-^; 
and jgaa are in units of 10'^°erg/s/h~^Mpc^, and jd may be larger than jgas due 
to the inaccuracy of the Schechter fit. 

We estimate that a is well constrained (±0.03) and that the product uqL* is 
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also fairly well constrained, but individual values (no, I/*) are poorly determined 
because L* is dependent on the quite uncertain highest luminosity clusters. To 
estimate the purely statistical uncertainty, we reanalyzed the z = data from 
the lower panel of Table la (total, integrated X-ray cluster properties) looking 
separately at two halves of the box. The fractional differences {\AQ/Q\) for 
Q = {a,L*,kT:r,no,jciJgas) were found to be (0.0062,0.90,0.62,0.74,1.31,0.76) 
respectively. The fact that even the integral j varies significantly between the 
two halves of the box reminds us again of the "cosmic variance". Our sample 
volume is not large enough to give us a robust estimate for the cosmic mean value 
of j. Note that the variances of all the quantities (except for a) in this model are 
larger than those in the CDM model (KCOR), presumably due to the fact that 
this CDM+A model has relatively more large-scale power than the CDM model, 
indicating that a larger box is clearly desirable. This, of course, is related to 
the virtues of open models: they typically have relatively more large-scale power 
and thus match better a variety of large-scale structure observational constraints 
(c/. , e.g., Efstathiou et al. 1993) than do = 1 models. 

We see that the cluster cores, as we have defined them, contain between 1/3 
and 1/2 of the total X-ray emission in the regions studied, comparable to the 
results of KCOR. The total cluster luminosity in the box is typically 3/4 of the X- 
ray emission from the box, the same level was found in the CDM model (KCOR). 
This fraction depends on energy, as one can see by comparing Tables lb and 4d. 
Since most of the very high temperature gas is in the central regions of clusters, 
the fraction of the total emission from identifiable clusters in the 2 — lOkeV 
region approaches unity. For the total luminosity the Schechter a parameter is 
approximately 1.55 (as compared to 1.50 in KCOR) with little evolution, and 
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CK, for the few keV bands is typically slightly flatter at 1.4. For o; < 2, as we 
have noted, most of the luminosity arises from bright clusters. Observationally 
a = 1.9 — 2.0; so this model provides a marginally better fit than does the 
standard Q = 1, CDM model. Allowance for line emission processes would 
increase the luminosity of the lower luminosity (low temperature) clusters more 
than the high luminosity (high temperature) clusters, still further improving the 
fit. Without a careful treatment of line emission, it would be premature to say if 
the fit to the slope of the luminosity function, as represented by a in this model, 
is adequate. 

The number density of bright clusters peaks at intermediate redshift, and 
the typical luminosity is (for small redshift) relatively constant, so there is a 
peak emissivity at approximately z = 0.2 — 0.5 for the several keV bands. Thus, 
crudely speaking, in this model one expects weak "positive" evolution until nearly 
z = 0.5 and then negative evolution thereafter. The peak occurs at significantly 
lower redshift in the open model than in the standard CDM model. Figure (5a) 
shows the comparison between the two models for the number density (per unit 
comoving volume) of clusters with rest frame luminosity > 10^'^erg/s. Figure 
(5b) shows the number density (per unit comoving volume) of clusters with rest 
frame luminosity (integrated over the entire frequency range ) > 10^3gj.g/s for five 
different models. We see that it is the overall normalization on the relevant scale 
(cTs) rather than Q which is the dominant factor on the evolutionary behavior of 
the bright cluster number density. Physically, the rapid evolution occurs at an 
epoch when most of the bright clusters are collapsing. After that, the evolution 
still goes on due to processes such as merging, but it is relatively mild. Comparing 
the COBE-normalized standard CDM and COBE-normalized CDM+A models, 
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negative evolution begins earlier and is stronger for the open model. Allowance 
for cosmological effects would make both models evolve more negatively, but these 
effects would be exaggerated in the CDM+A case due to the rapid increase in 
luminosity distance (with redshift) in cosmologies with a significant cosmological 
constant. The different evolutionary path of the cluster luminosity function is 
the most significant difference we have found between open and Q = 1 models. 
For example, Figure (5b) shows more than an order of magnitude decline in 
the number density of clusters with 0.5-4.5keV luminosity > lO^^erg/s between 
redshift zero and two for the open model, but it shows a slight increase for the 
standard CDM n = l model. 

We believe that the peak in emission seen at moderate redshift is real. The 
reasons for this are discussed in KCOR. The approximate Press-Schechter for- 
malism, which does not allow for a variety of effects, cannot easily mimic the full 
non-equilibrium hydrodynamic treatment. 

While the integrated X-ray emissivity evolves fairly slowly over the period 
surveyed in Table (la) with a flat maximum between z = 0.2 and z = 0.3, L* and 
no tend to evolve more rapidly and in opposite directions. The primary difference 
which we find between this model and the CDM model (KCOR) is the evolution 
of the bright end X-ray clusters. In the = 1 CDM model (KCOR) we find that 
L* is a monotonic decreasing function of redshift, while in this model we find 
L* peaks at 2; ~ 0.3. The reason for this difference is the dramatically different 
behaviors of the universal expansion at later times in the two models. We expect 
that in a lower O model, L* will clearly peak at an even higher redshift. The 
evolution of L* should be observable and provide a powerful discrimenant among 
cosmic models. 
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To highlight the negative evolution of the bright end of the luminosity func- 
tion, we listed in the fifth and sixth columns of Tables Ib-ld (sixth and seventh 
columns of Table la) the comoving density of clusters having luminosity greater 
than lO^^erg/s and lO^'^erg/s. For reasons stated earlier (based on our limited 
box size), we use the Schechter fit rather than direct counts to compute these 
columns. Comparing columns 4 and 6 (5 and 7 of Table la), we see that, although 
the total number density no of clusters increases with redshift (until ^ a; 1 — 2), 
the number density of the highest luminosity {L^ > lO^^erg/s) clusters decreases 
for z > 0.5. This is presumably one of the effects leading to the observational 
appearance of "negative evolution" . Statistical fluctuations in our results are still 
quite significant due to the limited box size, especially for n{Lx > lO'^'^erg/s). 

Redshift effects strongly exaggerate this tendency to observe negative evolu- 
tion, since higher redshift clusters tended to have lower temperatures (c/. column 
4 of Table la and Figures 6, 13), and both effects will reduce the energy observed 
by satellites measuring the X-ray flux in high energy bands. Note that the neg- 
ative evolution in the density of clusters with L > lO'^'^erg/s is more and more 
steep in Tables (a— > d) as one looks at higher energy bands, although statistical 
fluctuations are large. 

The emission- weighted temperature, T^, of each cluster is calculated and 
the distributions are shown in Figure (6). The arrow in each panel indicates the 
average cluster temperature (weighted by luminosity) at the given epoch. Also 
shown in Figure (6b) is the observed temperature function from Henry & Arnaud 
(1991) as the cross-shaded area. Note that the computed temperature function 
is somewhat lower than the observed one, but we think that the limited size of 
our simulation box caused an omission from the computation of the highest tem- 
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perature clusters; a bigger box is needed before we can have definite conclusions 
on this. We see that at all epochs the coolest clusters dominate the statistics 
(the turnover at low is presumably caused by our definition of minimun cell 
luminosity to constitute an X-ray bright cell), but the mean is determined by 
the high mass, high luminosity, high temperature end of the distribution. The 
mean temperatures, indicated by arrows, are included in column 3 of Table (la). 
We will return to the issue of temperature evolution later. Looking ahead to 
Figure (8a) , we see the strong correlation found between and total luminosity 
(clusters are shown at ^ = 0). 

Now let us turn again to the total cluster luminosity (vs. core luminosity) 
as shown in the lower panels of the tables and in Figures lb-4b, the quantity 
normally measured by satellite observations. The ratio of {jcl,tot/ jcl,core) is near 
2.0 ± 0.1 for the redshift range < 2; < 1; a similar value was found in KCOR. 

We can also roughly estimate the effective radii of the clusters by assuming 
that the emission has a profile 

i = ^2 (2) 

and determining, from the ratio of the luminosity (integrated over frequency) 
of the central cell to the total cluster luminosity, the value of Vx which would 
produce this ratio. We show in Figure (7) the radii determined in this fashion. 
The peaks seen in the panels of Figure (7) are, of course, artificially induced 
by our cell size of 0.31/i~^Mpc, but the distribution to larger radii should be 
reasonably accurate. Arrows indicate the luminosity-weighted average values. 
Since brighter clusters tend to be resolved, these numbers should be reliable, but 
the arrows are uncomfortably close to the peaks of the curves, indicating that 
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the most luminous clusters may be somewhat unresolved. We see a weak trend 
of increasing size with increasing time, which is in the theoretically anticipated 
direction. Longer wavelengths became nonlinear later, producing larger clusters, 
and smaller clusters merge to produce larger clusters with increasing time. 

Now in Figures (8a,b) we show the scatter plots of (Tx,ra;) vs. (inte- 
grated over frequency). We see that there is a clear correlation between and 
Ta;. But we do not see any strong correlation between and r^- The best 
fit lines (dashed) indicate a slope of (0.36 ± 0.01,0.32 ± 0.01) for (core, total) 
cluster region. In KCOR we found (0.394 ± 0.001, 0.375 ± 0.001) for the CDM 
model. The observed correlation [cf. Figure (8a) shown as cross-shaded area] be- 
tween Tx and total cluster luminosity is log^g Tx{keV) — log^g (4.2^0^8) + (0.265 ± 
0.035) logiQ(/i^L44), according to Henry & Arnaud (1991). In the region where 
a comparison can be made, the agreement with observations is good. 

Next, we address temperature variations within clusters. Given our lim- 
ited resolution, there is little that can be accurately determined on this issue 
from our simulations, but we are able to compare the central cell (Volume= 
3.1 X 10~^/i~^Mpc'^) with the surrounding cells (3^ — 1^) (volume of size 5.9 x 
10~-^/j.~^Mpc^) and the cells surrounding these cells (5^ - 3^, vol= 3.0/j.~^Mpc^). 
We define the ratio of the inner cell to the next cube as Tc/Tgheii (volume- 
weighted average) and show the scatter diagram in Figure (9). No trend is seen 
with luminosity, and the median value, indicated by the dashed line, is 1.4. The 
cluster gas deviates significantly from isothermality, with a 5%— 10% temperature 
decline typically found by a radius of 0.4 — 0.5/j.~^Mpc but a sharp fall-off is in- 
dicated for radii l/i~^Mpc. In Figure (10) we compare the (luminosity-weighted) 
temperatures found in the three regions noted above, normalized to the tem- 
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perature in the central cell. The large dispersion is indicated by the error bar. 
Figures (9) and (10) show no significant differences from the analogous figures in 
the = 1 case (KCOR). 

We show in Fig. 11 and 12 the evolution of cluster core radius (in metric, not 
comoving, units) and temperature (both luminosity-weighted) as open circles. 
Solid dots are from the KCOR Q — 1 computation. Also shown are the fits 
analytically predicted by Kaiser (1986) i?a; oc (1 + z)~'^ and T^; oc (1 + z)~^ for 
n = 1, CDM model (solid line) from KCOR and this CDM+A model (dotted 
line) (note that Kaiser's prediction is valid for = 1 models only, but is shown 
also for the CDM+A model for comparison purposes). Examination of these 
figures indicates the radius changes in a way qualitatively similar to that found 
for standard CDM with an ~ 40% decline in radius from redshift zero to a look- 
back of redshift z — 0.5. But there is an overall difference in that the radii are 
smaller in the CDM+A model. The difference, while small (—23%), is significant, 
especially since the lower bound produced by our finite numerical resolution is 
less than a factor of two below the mean value. Since real, observed clusters are 
smaller than the computed clusters, this change is also in the desired direction. 
The differences between the two models are dramatic with regard to temperature. 
First, as noted, the mean temperature is a factor of 3.5 lower in the open model 
and the evolution is far slower (—33% to ^ = 0.5 rather than —48% for = 1). 

Thus we see that the differences between the two models at low redshift 
are significant and very important, due to freeze-out at low redshift in open 
models. This temperature evolution should be detectable even with a relatively 
"soft" X-ray instrument such as ROSAT. It will be able to discriminate between 
f2 = 1 models and Q < 1 models. Furthermore, detailed and direct numerical 
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simulations combined with more observations might provide a way to constrain 
fl, or A, or combinations of Q and A. 

Figures (13) and (14) show the evolution in the background radiation field 
in two additional ways. The first shows what a comoving observer would have 
measured at various redshifts. But below IkeV the results are unreliable because 
of both the omission of line emission and the omission of IGM absorption. The 
second, Figure (14), shows the fractional contribution to the background seen by 
an observer at 2; = in several bands that were produced at various epochs (in 
integral form). The important point to note is that most of the X-ray background 
(especially in the harder bands) that we see locally were produced at relatively 
recent {z < 0.5) epochs. This is a consequence of many things, prominant among 
them are the redshift factors that dilute the observable effects of emission at high 
redshift. But there is a major difference between the two models. In this open 
model, one-half of the 2 — lOkeV brightness of the sky comes from redshift less 
than z — 0.3, whereas, in the Q — 1 case (KCOR), the median point is much 
closer, at redshift z = 0.2. 

Finally, let us take a slightly different route to address the issue of bias of 
gas relative to the mass: does the gas in dense regions, like clusters of galaxies, 
fairly represent the underlying mass, or is it "biased" or anti-biased? This is 
a question with great cosmological significance. If we know the ratio of gas (+ 
galaxies) to total matter in the clusters by direct observations, and we know, from 
light element nucleosynthesis the global baryon density, then we can divide the 
second number by the first to obtain the global matter density and to compare 
with the cherished critical density. This old line of argument has been carefully 
re-examined recently by White (1992) and also reanalyzed by Babul & Katz 
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(1993) and others. 

The argument depends on knowing whether or not Ph/ Ptot varies significantly 
from place to place and, in particular, whether this quantity will be found near its 
average value in the high density regions, where it can most easily be measured. 

Our possibly counter-intuitive results are shown in Figure (15), where we 
plot the ratio {pigm / Ptot) vs. {ptot / {Ptot)) smoothed by a gaussian of radius 
l/j.~^Mpc. At any value for the total density, there is a wide range of possible 
values of pioM, but the high density regions actually have a lower than average 
ratio of baryons to total mass. We see that, in the high density clusters, where 
Ptot/ < Ptot > approaches 10^, the gas is under represented by a factor of 1.1. 
(In KCOR we found a factor of 1.7 for the CDM model.) Combining this factor 
with its global mean, we obtain that (^^)c/ = 0.083 ±0.007, which is consistent 
with observations [{^^ff)obs = 0.097 ± 0.019 for the Coma cluster for the given 
value oih, cf. White 1991]. 

We can use the comparison with KCOR to address the question of whether or 
not the small anti-bias found is real or due merely to numerical errors spreading 
out the gas density more than it does the dark matter density. In this simulation 
the core radii of the clusters are smaller on average than in the f2 = 1 case, 
putting them in better agreement with observations but, as noted, closer to our 
grid spacing. Thus numerical diffusion in this case should be a larger (relative) 
effect here than in KCOR, and the anti-bias, if it were due to numerical diflFusion, 
would be larger. But it is smaller, indicating that this is not the case. Our best 
guess is that the cause of the anti-bias is related to outward propagating shocks 
in the transient formation phase. Since such time-dependent effects will be less 
in open models, the anti-bias should be less here - as it is. 
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Now, we treat the same problem in a somewhat different way concentrating 
on the identified clusters. For each simulated X-ray cluster, we compute the gas 
mass and total mass within a sphere of radius l/j.~-'^Mpc (centered at the X-ray 
cluster center) and in Figure (16) we show the ratio (open circles) of the two 
masses within this sphere as a function of density of the sphere relative to the 
mean. The dotted line is the best log-log straight line fit for the open circles, and 
the solid line is the fit weighted by the luminosity of each cluster. Also shown, 
as the dashed line, is the global mean of the ratio. Presenting the same infor- 
mation in an alternative way, we show in Figure (17) the histogram of the ratio, 
both number- weighted (thin solid histogram) and luminosity-weighted (thin dot- 
ted histogram) in the CDM-I-A model. The heavy solid histogram (arbitrarily 
normalized to have a similar peak height) indicating the observational situation 
is adapted from Jones &: Forman (1992). We see that there is a trend that poor 
clusters are relatively gas poor, which is consistent with observations (c/. , e.g., 
Jones & Forman 1992). The luminosity- weighted fit is closer to (but still below) 
the global mean with an anti-bias in the range 0.85 — 0.92. We see that the 
median of the computed ratio is in agreement with the median of the observed 
ratio, whereas in the CDM Q = 1 model the computed ratio is lower by a factor 
of 2-3. Improving the observations will probably narrow the heavy histogram, 
and increasing the dynamical range of our simulation box will widen the thin 
histograms. Both expected improvements should make the agreement better. 

4. CONCLUSIONS 

In Figures (1-5) we show the evolution of the cluster luminosity function ex- 
pected in the CDM-I-A model. In the range of parameters where there is greatest 
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overlap between observed and computed quantities (0 < 2; < 1, 10^°erg/s< L^^ < 
lO^^erg/s), little evolution is seen (for comoving observers) in any of the com- 
puted bands aside from a decline in the number of brightest sources (somewhat 
uncertain due to our limited box size) and a modest increase (by about a factor 
of two) in the luminosity function for fainter objects. A similar behavior was 
also found in the CDM model (KCOR). But the likely explanation for this is 
different from that in the CDM case. While in the standard CDM model, it is 
largely coincidental and due to the balancing of two effects: new breaking waves 
increase the luminosity density but mergers decrease it. In this CDM+A model it 
is primarily due to the late time freeze-out of formed clusters. The effect is seen 
most clearly in Figure (5), where we look only at the evolution of the brightest 
clusters, and see an earlier decline (with increasing redshift) of the open model 
than the = 1 models. There is an order of magnitude decline in the number 
density of clusters having > lO^^erg/s (in the harder energy bands) in the 
redshift interval 2; = ^ 2 for this model but an increase in the standard = 1 
CDM model. The strong negative evolution found in this paper in an open model 
would of course be greatly enhanced for observers (using ftxed bands and inten- 
sity limits) by cosmological effects which are especially strong in models with a 
cosmological constant (Carroll, Press, & Turner 1992) 

Figures 6 and 7 show rates of change in other quantities, the temperatures 
and clusters radii. These important trends are summarized in Figures 11 and 12 
where we see a factor of 2-3 decline in both these quantities by redshift 1. But 
more interesting is the nearly constant mean temperature (luminosity-weighted) 
in the redshift range z = ^ 0.3 in the CDM-I-A while in the CDM model we 
see a sharp decrease in the temperature from z = 0.0 to z = 0.3. Equivalently 
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put, L* of Schechter fits to the computed luminosity functions peaks near z = 0.3 
in the CDM+A model but it increases monotonically until z — in the CDM 
model. Both trends should be detectable even with a relatively "soft" X-ray 
instrument such as ROSAT. They might provide powerful tests for = 1 models 
and < 1 models. Also the actual values of the cluster temperature are much 
lower in the open model and in better agreement with observations. 

This CDM+A model provides an adequate fit to the observed bright X-ray 
cluster luminosity function and X-ray background; however, we found that there 
would be too many bright X-ray clusters produced and too much integrated 
background X-ray intensity in the COBE-normalized standard CDM model. 

We find a slight anti-bias (~ 10%) of gas relative to the mass in dense regions 
like the clusters of galaxies, but the model is self-consistent in the sense that the 
computed (g^)c/(= 0.083 ± 0.007) is consistent with observations [(gff)o6s = 
0.097 ± 0.019 for the Coma cluster for the given value of h, cf. White 1991], 
whereas there was a gross inconsistency in the VL = 1 case. 

In sum, the model is significantly different from that obtained from the stan- 
dard ri = 1, CDM simulation, and, with regard to all measurable quantities that 
we have compared to observations [N{Lx)., (Rx), (Tx), dlnTx/dlnL^, Pgas/Ptot\i 
it is not only a better fit to observations than standard il — 1^ CDM model, but 
also is an adequate representation of them. The predicted evolutionary differ- 
ences between open and Vt = 1 models are sufficiently great to allow definitive 
tests by current or planned X-ray satellite observations. 
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FIGURE CAPTION 



Fig. 1- Figure (la): the X-ray cluster bremsstrahlung luminosity (from central < 0.5h~^Mp(^ 
regions) function integrated over the whole frequency range at five different red- 
shifts z = (0,0.2,0.5,0.7,1.0). Figure (lb): the X-ray cluster bremsstrahlung 
luminosity (from < 1.0/i~^Mpc region) function integrated over the whole fre- 
quency range (filled dots) at the same five different redshifts. The cross-shaded 
area shows the observations (Henry & Arnaud 1991, {3.lt.tl x lO~^h^Mpc~^h'^ [i^44(bol)] "^j x| 
[/i2L44(bol)]-i-S^±0-4). 

Fig. 2~ Same as Figure (1) but for the luminosities integrated over 0.3 — 3.5keV frequency 
bin. 

Fig. 3- Same as Figure (1) but for the luminosities integrated over 0.5— 4.5keV frequency 
bin. 

Fig. 4- Same as Figure (1) but for the luminosities integrated over 2 — lOkeV frequency 
bin. The cross-shaded area in (4b) indicates observations (Henry 1992). 

Fig. 5- Figure (5a) shows the comparison between the two models for the number density 
(per unit comoving volume) of clusters with rest frame luminosity > lO^^erg/s 
for three X-ray bands. Figure (5b) shows the number density (per unit comov- 
ing volume) of clusters with rest frame luminosity (integrated over the entire 
frequency range) > 10^"^erg/s for five different models. 

Fig. 6- The X-ray cluster temperature {T^, emission- weighted temperature) function 
at six different redshifts z = (0,0.5,1,2,3,5). Arrows indicate the luminosity- 
weighted average temperature Tx at each epoch. In the first {z = 0) panel in 
(5b) the cross-shaded area is the observed temperature function from Henry & 
Arnaud (1991) [(1.81°;^ x 10-^h^Mpc-^keV-^){kT)-'^-'^'^'^-% 



Fig. 7- The X-ray cluster effective radius (r^) distribution [cf. equation (2)]. Arrows 
indicate the luminosity-weighted effective radius at each epoch. 

Fig. 8- Figure (8a) shows the scatter plot of Tx vs z = 0. The cross-shaded area 

line indicates the observations of Henry & Arnaud (1991). The dashed line is the 
best fit of the simulation results. Figure (8b) shows the scatter plot of Vx vs 
dX z — 

Fig. 9- The ratio of the central cell temperature to the temperature of its surrounding 
shell (~ one cell thick) as a function of Ltot- 

Fig. 10- We compare the (luminosity-weighted) temperatures found in the three regions 
(central cell, the shell surrounding the central and the next outer shell) and 
normalized to the temperature in the central cell. Departure from isothermality 
increase significantly for hr > O.SMpc. Note the errorbars are la variance. 

Fig. 11- The average cluster core radii in physical units as a function of redshift for 
clusters with luminosity in the 0.5 — 4.5keV band greater than lO^^erg/s for 
CDM-I-A model (open circles, this paper) and the standard = 1, CDM model 
(solid dots, KCOR). The best fit evolutions of the form oc {l + z)~'^ are shown 
as a solid curve for the CDM model and a dotted curve for the CDM-I-A model. 

Fig. 12- The average cluster temperature as a function of redshift for clusters with lumi- 
nosity in the 0.5 — 4.5keV band greater than 10^"^erg/s for CDM-I-A model (open 
circles, this paper) and the standard fl = 1, CDM model (solid dots, KCOR). 
The best fit evolutions of the form Tj; oc (1 -|- z)~^ are shown as a solid curve 
for the CDM model and a dotted curve for the CDM-I-A model. Temperatures 
in the CDM-I-A model are lower and evolution less than in the standard 0, = 1, 
CDM model, and tend to freeze out at lower redshift while we see a dramatic 
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increase in the CDM model approaching z = 0. 

Fig. 13- The mean radiation intensity at six epochs, z = 5 (solid line), z = S (dotted 
line), z = 2 (short, dashed line), z = 1 (long, dashed line) z = 0.5 (dotted-short- 
dashed line) and z = (dotted- long-dashed line). The box in the middle shows 
the obseravtional data by Wu et al. (1991). Neither line absorption nor emission 
has been allowed for in this figure. 

Fig. 14- The distribution functions of four presently observed X-ray bands as a function 
of redshift (in integral form). 

Fig. 15- The ratio Pigm/ Ptot as a function of ptot/ < Ptot >• Results are smoothed by 
a gaussian window of radius lh~^Mpc. The global mean value of Pigm/ Ptot is 
shown by the dashed line. Note that in the highest density regions the gas is 
under-represented, "anti- biased" , by a factor of about 1.1 (which is less than 1.7, 
found for n = 1 CDM model in KCOR). 

Fig. 16- The ratio p gas/ Ptot as a function of ptot/ < Ptot > within a radius of lh~^Mpc for 
each identified cluster (open circles). The dotted line is the best log-log straight 
line fit for the open circles and the solid line the fit weighted by the luminosity 
of each cluster. Also shown as the dashed line is the global mean of the ratio. 
We see that there is a trend that poor clusters are relatively gas poor. 

Fig. 17- The histogram of the ratio both number- weighted (thin solid histogram) and 
luminosity-weighted (thin dotted histogram) in the CDM-I-A model. The heavy 
solid histogram indicating the observational situation is adapted from Jones & 
Forman (1992). We see that there is a trend that poor clusters are relatively 
gas poor, which is consistent with observations (c/. , e.g., Jones & Forman 1992). 
The luminosity-weighted fit is closer to (but still below) the global mean, with an 
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anti-bias in the range 0.85 — 0.92. We see that the median of the computed ratio 
is in agreement with the median of the observed ratio, whereas in the CDM fl = 1 
model the computed ratio is lower by a factor of 2-3. Improving the observations 
will probably narrow the heavy histogram, and increasing the dynamical range of 
our simulation box will widen the thin histograms. Both expected improvements 
should make the agreement better. 
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Table la. Parameters of Schechter fits for the X-ray cluster luminosity function 

integrated over the entire frequency range. 
X-ray Cluster Core Luminosity (< Q.5/t~^Mpc) 



kJC 








■n(L-> 10^3\ 


n(Tj > 10^^") 


jcl 




n '^4 


1 7Q 


4 59 




7 4x1 n~2 


n 045 


1 fin 


1 16 


1 51 


9 3Q 


4 7 


94 


061 




n 97 

VJ. Z( 1 


1 05 


1 n 7 


4 9 


1 1 y in~2 




1 54 


15 


84 


19 1 


3 1 


4.8 X 10"^ 


055 


1.53 


0.11 


0.64 


27.7 


2.6 


4.1 X 10"^ 


0.057 


1.58 


0.043 


0.31 


37.9 


0.28 


8.8 X 10-12 


0.034 


1.93 


0.22 


0.14 


1.14 


0.38 


2.1 X 10"^ 


0.035 


1.90 


0.11 


0.054 


0.37 


3.1 X 10-2 


2.3 X 10-^ 


0.0039 


X-ray Cluster Total Luminosity {<\h -"^Mpc) 


OL 


LI (10^4) 


k% (keV) 


no 


n(L > 10^3) 


n(L > 10^4) 


jcl 


1.56 


1.74 


1.55 


2.56 


6.6 


0.54 


0.090 


1.56 


1.86 


1.29 


3.14 


8.5 


0.73 


0.12 


1.54 


0.59 


0.91 


9.64 


9.2 


0.20 


0.11 


1.50 


0.29 


0.72 


20.4 


8.6 


3.2 X 10-2 


0.10 


1.48 


0.15 


0.58 


36.3 


6.0 


1.0 X 10"^ 


0.093 


1.59 


0.082 


0.30 


30.4 


1.5 


1.1 X 10-^ 


0.054 


1.78 


0.046 


0.15 


15.1 


0.11 


9.4 X 10-12 


0.029 


1.90 


0.020 


0.020 


3.56 


3.7 X 10-^ 


9.8 X 10-24 


0.0068 



Here L* and L are in units of lO^^erg/s; no, n(L > 10^^) and n(L > lO'^^) are in 
units of 10-^/i^Mpc-^; jd and jgas are in units of lO'^^erg/s/h-'^Mpc'^, and jd 
may be larger than jgas due to the inaccuracy of the Schechter fit. 



Table lb. Parameters of Schechter fits for the X-ray cluster luminosity function 

in 0.3 — 3.5keV band 











nCL > 10^3) 


n(L> IQ^'^) 


Jcl 


J gas 


J 


1 45 


0.22 




1 7 

-Lit 


2.4 X 10~^ 


021 


060 


2 


1 48 


25 


6 05 


2.1 


4 7 X 10"^ 


026 


070 




1 45 


26 


q 32 


2.4 


9 3x1 0~^ 


039 


075 


7 


1.46 


0.24 


9.66 


3.2 


6.3 X 10-^ 


0.038 


0.069 


L 


1.35 


0.22 


8.40 


2.4 


4.0 X 10-3 


0.026 


0.059 


I 


1.38 


0.14 


5.25 


0.78 


1.0 X 10-^ 


0.011 


0.020 


I 


1.60 


0.13 


3.02 


0.38 


1.9 X 10"^ 


0.0087 


0.003^ 


5 


1.90 


0.010 


0.03 


6.4 X 10-9 


4.4 X 10-26 


2.8 X 10-^ 


2.6 X 10 



X-ray Cluster Total Luminosity {<lh -"^Mpc) 





a 


T * 


no 


n(L > 10^3) 


n(L > 10^4) 


jcl 


3 gas 


) 


1.42 


0.44 


6.37 


4.1 


6.0 X 10-2 


0.043 


0.060 


2 


1.43 


0.48 


6.94 


4.9 


8.5 X 10-2 


0.052 


0.070 


5 


1.37 


0.32 


14.2 


6.3 


4.2 X 10-2 


0.065 


0.075 


7 


1.31 


0.13 


30.5 


4.1 


3.7 X 10-^ 


0.052 


0.069 


L 


1.36 


0.16 


22.4 


4.2 


1.3 X 10-3 


0.050 


0.059 


I 


1.24 


0.13 


9.59 


1.3 


1.3 X 10"^ 


0.015 


0.020 


? 


1.60 


0.35 


0.77 


2.5 X 10-3 


6.0 X 10-^6 


0.0060 


0.003^ 


5 


1.90 


0.01 


0.050 


1.1 X 10-8 


7.4 X 10-26 


4.8 X 10"^ 


2.6 X 10 
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Table Ic. Parameters of Schechter fits for the X-ray cluster luminosity function 

in 0.5 — 4.5keV band 









Tin 


niL> 10^3) 


n(L> IQ^'^) 


Jcl 


J gas 


J 




0.21 




1 7 

-Lit 


2.2 X 10~^ 


019 


053 


2 


1.44 


n 26 


5 74 


2.1 


5 8x1 0"^ 


024 


059 




1 46 


57 


4 65 


4 


9.4 X 10~2 


044 


061 


7 


1.39 


0.16 


12.1 


2.2 


6.7 X 10-^ 


0.028 


0.053 


L 


1.34 


0.13 


10.9 


1.4 


1.2 X 10-4 


0.019 


0.042 


I 


1.33 


0.054 


9.30 


0.18 


7.2 X 10-^° 


0.0068 


0.012 


I 


1.58 


0.025 


5.2 


3.5 X 10"^ 


2.9 X 10-20 


2.7 X 10-4 


1.5 X 10 


5 


1.65 


0.09 


1.60 


9.2 X 10-2 


1.7 X 10-^ 


0.004 


5.4 X 10 



X-ray Cluster Total Luminosity {<lh -"^Mpc) 





a 


T * 


no 


n(L > 10^3) 


n(L > 10^4) 


jcl 


3 gas 


) 


1.36 


0.44 


6.46 


4.0 


6.4 X 10-2 


0.040 


0.053 


2 


1.42 


0.58 


5.23 


4.4 


0.11 


0.047 


0.059 


5 


1.40 


0.52 


7.92 


5.9 


0.13 


0.061 


0.061 


7 


1.32 


0.14 


20.6 


3.1 


4.6 X 10-4 


0.038 


0.053 


L 


1.31 


0.12 


22.4 


2.6 


1.3 X 10-4 


0.035 


0.042 


I 


1.38 


0.051 


11.7 


0.18 


2.4 X 10-^° 


0.0086 


0.012 


? 


1.58 


0.035 


5.2 


1.7 X 10-2 


4.3 X 10-^^ 


3.8 X 10-3 


1.5 X 10 


5 


1.43 


5.95 


0.47 


3.1 X 10-2 


8.8 X 10-s 


0.044 


5.4 X 10 
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Table Id. Parameters of Schechter fits for the X-ray cluster luminosity function 

in 2 - lOkeV band 



X-ray Cluster Core Luminosi 


ty (< 0.5/i-iMpc) 




Uc 




IIQ 




nfr, ^ 1 n^^'i 


Jcl 


jgas 




1 99 


n nss 

U.UOO 




U.U 1 


9 A V 1 


n ni 1 

U.Ul 1 




9 


1 "^7 


n 1 fi 




n Q9 


9 Q V 1 


n ni 1 


n 09/1 


K 

o 




n 1 

\J.i-0 


U. 1 u 


1 1 


9 ^ V 1 


n nizL 


n ni Q 


7 




17 


3 1 3 




q N/ in~^ 


0074 


01 3 




1.20 


0.10 


3.88 


0.34 


4.4 X 10-6 


0.0045 


0.0075 




1.17 


0.081 


0.79 


4.7 X 10-2 


7.4 X 10-s 


7.0 X 10-^ 


7.8 X 10 




1.13 


0.018 


0.14 


2.9 X 10-5 


5.1 X 10-23 


2.8 X 10-5 


3.1 X 10 




1.80 


0.004 


0.01 


1.8 X 10-16 


1.3 X 10-26 


1.8 X 10-6 


2.5 X 10" 


X-ray Cluster Total Luminosity {< Ih -"^Mpc) 




a 


T * 


no 


n(L > 10^3) 


n(L > 10^4) 


jcl 


Jgas 




1.19 


0.18 


9.40 


2.0 


1.7 X 10-3 


0.020 


0.026 


2 


1.35 


0.30 


5.88 


2.4 


1.3 X 10-2 


0.024 


0.024 


5 


1.31 


0.25 


5.01 


1.7 


5.1 X 10-3 


0.016 


0.019 


7 


1.32 


0.20 


2.82 


0.71 


8.1 X 10-^ 


0.010 


0.013 




1.13 


0.17 


3.14 


0.63 


4.4 X 10-^ 


0.0058 


0.0075 




1.10 


0.10 


1.25 


0.11 


1.8 X 10-6 


1.3 X 10-3 


7.8 X 10 




1.03 


0.036 


0.14 


1.0 X 10-3 


1.7 X 10-15 


5.2 X 10-5 


3.1 X 10 




1.80 


0.005 


0.01 


3.8 X 10-1^ 


1.5 X 10-26 


2.3 X 10-6 


2.5 X 10" 
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